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A minimal model is exactly solved for electron spin transport on a helix. Electron transport is 
assnmed to be supported by well oriented pz type orbitals on base molecules forming a staircase of 
definite chirality. In a tight binding interpretation, the SOC opens up an effective ttz — rz coupling 
via interbase pz:,y —Pz hopping, introducing spin coupled transport. The resulting continuum model 
spectrum shows two Kramers doublet transport channels with a gap proportional to the SOC. Each 
doubly degenerate channel satisfies time reversal symmetry, nevertheless, a bias chooses a transport 
direction and thus selects for spin orientation. The model predicts which spin orientation is selected 
depending on chirality and bias, changes in spin preference as a function of input Fermi level and 
scattering suppression protected by the SO gap. We compute the spin current with a definite helicity 
and find it to be proportional to the torsion of the chiral structure and the non-adiabatic Aharonov- 
Anandan phase. To describe room temperature transport we assume that the total transmission is 
the result of a product of coherent steps limited by the coherence length. 


I. INTRODUCTION 

Spin selectivity is an exciting prospect for many tech¬ 
nological applications that require spin imbalance as a re¬ 
source for information storage and processing [I|, Q. Spin 
polarized currents, are traditionally derived from opti¬ 
cal orientation^ when material-optic interphases are re¬ 
quired or ferromagnetsQ involving gross power sources 
to control. Spin polarized currents derived by more sub¬ 
tle means are then desirable in the spin-torque driven 
writing on magnetic surfaces, valve devices and more gen¬ 
erally in spin based information processing gates. 

There is now well established experimental evidence 
that chiral molecules have the tantalizing property of 
polarizing electrons scattered either through 2D films or 
driven through individual chiral molecules These 

represent two limiting cases for transport: For 2D films, 
electrons are emitted from a metallic surface and are 
asymptotically free, having and energy above the film po¬ 
tential barrier. On the other hand, for transport through 
a single molecule, electrons tunnel between metallic con¬ 
tacts driven by a potential bias or a two terminal set 

up0,ll. 

The former, scattering experimental setup, is un¬ 
derstood theoretically in terms of multiple scattering 
through the chiral target in the presence spin-orbit cou¬ 
pling (SOC)[l-[ll|. The theoretical model requires inter¬ 
ference between scalar potential events and spin-coupled 
scattering. Chiral structure and SOC define an en¬ 
ergy window in which the electron spin is effectively 
polarized[l3, [13 ™ tti® range of 1-10 eV. A higher en¬ 
ergy limit is required for multiple coherent scattering and 
electron probing of the chiral structure, while a minimal 
energy is required for the SO coupling to be significant. 
The spin-orbit coupling strengths required are compati¬ 
ble with those expected form the atomic carbon cores, in 
the energy scale of the meV. Room temperature polari¬ 


sation is proposed to be due to incoherent superposition 
of short coherent length nolarizationfl^ ■ 

On the other hand in the experiments performed W 
biasing a DNA molecule attached to metallic contacts^ 
have been addressed theoretically by approaches based 
on quantum transport [3 flsI - fT^ in an chiral electric field 
generated by the charge distribution of the molecule, or 
helical potentials. The electrons are considered to follow 
paths on the helix sampling the chiral electric field. An 
important feature of the model is the relation between the 
field intensities producing SOC and the electron mobility. 
The models predict a strong polarisation effect in spite of 
the small SOC because the low mobility enhances the res¬ 
ident times of electrons in the potential. Considerations 
of the limited coherence length due to room temperature 
operations have been discussed in refill. 

We aim here to theoretically describe the single 
molecule experiments in more detail: For electrons with 
energies below the molecular barrier, a current is driven 
by a potential difference, through single DNA molecules 
in an STM set upQ- The experimental results consist 
of a spin dependent barrier to tunnelling of electrons 
through the chiral molecule. A preferred spin direction 
results in a larger current than the opposite spin direc¬ 
tion. There is a gap between barrier heights that is inde¬ 
pendent of the molecule length and is thus proportional 
to the SO coupling. A striking characteristic of the TV 
curve is that it is anti-symmetric {V —>■ —V), and the 
preferred spin filtering direction in one bias direction is 
opposite to the preferred spin of the other Q. Finally, the 
bias needed to produce the same spin polarised current 
increases with molecule length. 

Here we present a model for the tunnelling of electrons 
through a chiral molecule which can be solved exactly 
explaining many experimental features. The model in¬ 
cludes the intrinsic SOC whose source are the atomic 
cores of the carbon atoms sitting on the bases of the 
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DNA structure (with a definite orientation in space). The 
strength of such a coupling for tunnelling electrons can 
be readily derived from a tight binding approach, and in¬ 
volves both the intrinsic atomic value and contributions 
from the coupling between nearest neighbour bases sites 
of the DNA molecule. We can anticipate the results for 
this model by performing a symmetry analysis: Both the 
SO coupling and the helix break inversion symmetry, but 
time reversal symmetry is preserved. This implies that 
the spectrum of the spin-orbit active helix should be com¬ 
posed of Kramers degenerate doublets, separated by the 
effective spin-orbit coupling gap. The quantum numbers 
of the helix of definite chirality comprise the kinetic en¬ 
ergy index, the rotation sense of the electron and its spin. 
Each Kramers doublet preserves time reversal symmetry, 
so that they comprise both rotation and both spin quan¬ 
tum numbers. On choosing a bias direction, only the 
channel (one of the two states in the doublet) for that 
sense of propagation is selected and has an associated 
spin. Thus, the bias breaks time reversal symmetry by 
only populating one channel and spin selectivity results. 


This paper is organised as follows: We first derive the 
Hamiltonian for a continuous one dimensional helix of a 
fixed number of turns and chirality, in the presence of 
spin-orbit coupling, whose source is a local atomic core 
electric field in the 2 direction. Electrons are constrained 
to follow the helical path on the corresponding eigenchan- 
nels. Two energies are defined; that of the free electron 
problem and the SO energy whose ratio determines the 
adiabaticity of spin transport (l^. Following, we obtain 
the channel energies and the corresponding exact eigen¬ 
functions as a function of two quantum numbers (current 
direction and spin) and a chirality index. We then show 
that the spectrum implies that there are always two dou¬ 
bly degenerate levels and each degenerate pair combines 
time reversed states. Whenever one biased direction is 
chosen, time reversal is broken and a preferred spin is 
filtered for that direction. The opposite spin state in the 
same direction is at a higher energy so that it requires 
a higher bias to be occupied, a feature that suppresses 
backscatteringQ. The resulting spin selected transport 
is evidenced through the computation of the spin current 
along the helix, leading to spin accumulation observed in 
the experiments. 


The hnal section addresses the problem of the SOC 
strength, where we report a tight-binding result which 
might address the strength of the spin-orbit energy gap 
observed. We also argue that the tunnelling through the 
DNA molecule occurs through an incoherent sequence di¬ 
viding the molecule into coherence-length long segments 
to explain the scaling of the barrier potentials. We end 
with the conclusions. 


II. DERIVATION OF THE INTRINSIC 
ELECTRIC FIELD HAMILTONIAN ON A HELIX 


The model Hamiltonian we propose consists of an elec¬ 
tron confined on a helix of N turns of radius a and pitch 
b as seen in FiglT] An internal electric field, of atomic 
origin, is assumed to exist as the source of the electron 
spin-orbit coupling. In an atomic model for the helix, the 
SO interaction offers a Pz to Px,y hopping route, first or¬ 
der in the atomic intrinsic SO coupling and thus, at least 
in the meV energy range. Assuming a motion strictly on 
the helix through a sequence of nearest neighbour states, 
we have the following Hamiltonian (see also ref. [ 13 ) 

H — H^inetic H" H30 

= +^'^) + {a/h){px(Ty-py(Tx) ( 1 ) 

where the first term is pure kinetic energy, while the sec¬ 
ond term is the SO coupling for an intrinsic electric field 
in the 2 direction. The strength of the electric field is em¬ 
bedded in the parameter a = eb?E/(Am?c^) which has 
dimensions of energy times a length. Nevertheless, this 
electric field is a nontrivial quantity to assess, since it not 
only contains a measure of the field felt by the electrons in 
their excursion to their nuclei, but also a quasi-resonant 
coupling contribution to the neighbouring statesfT^ [l^. 
We will discuss these contributions later. 

Due to the symmetry of the problem it is preferable to 
work in cylindrical coordinates, and after putting care¬ 
ful attention to hermiticity issues, well discussed in the 
literature [ 13,113 , we arrive at 

^ pfp _ I _ 

^ 2m*(a2 -t 62) “ h{a? + ^ 2(a? + b'^)'^ 

( 2 ) 

where Py, = m*(o? -I- b'^)^ = —ihdy,, a is the helix ra¬ 
dius, b is the helix pitch, and cTp = cr^ cos ip + (Ty sin ip 
and (Tip = —cr 2 ; sin -|- (Ty cos ip. The helix curvature is 
given by the ratio k = a/(a^ -|- 6^) and the torsion is 
T = bl((j? + b‘^). The momentum in the z direction will 
then be Pz = rp^. m* is the electron effective mass, 
which assumes the carbon states form a narrow band of 
states. It is physically convenient [T3| to identify two dis¬ 
tinct frequencies ujq = hlm*{a? + b"^) related to the free 
electron kinetic energy and tuso = 2aa/h(a?+b‘^) propor¬ 
tional to the helix curvature and the SO coupling. One 
can then simplify the the Hamiltonian in Eql2] into the 
simple quadratic form 




( 3 ) 


The ansatz for the wave function for spinfull electrons is 
on a helix of N turns with hard wall boundary conditions 
at the helix ends 


Xj/A.C _ ^iX{n/2N)<p 


/^A,Ce-iv/2\ 


( 4 ) 
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FIG. 1. Energies of electrons on a counter-clockwise helix 
(C = -t-1) as a function of the SO interaction strength cuso- 
The energy is in units of hcuo. Each level is doubly degenerate 
with the index labels indicated. When cuso = 0 each energy 
is fourfold degenerate as expected since the time inversion 
symmetry turns into independent space and spin inversion 
symmetries. A indicates the energy gap between Kramers 
degenerate states. The inset shows the SO magnetic field and 
the orientation of the spin for an eigenstate of the spin-orbit 
coupled system. 


where A = -|-1(—1) labels for the counter clockwise, 
(clockwise) electrons, s = ±1 labels the spin and C = ±1 
labels the chirality of the helix. Although the chiral in¬ 
dex does not appear in the Hamiltonian, it chooses the 2 
direction of propagation for a particular A index. Using 
the wave function ansatz one can derive the exact energy 
of the model as 


n,s 



(5) 


where n is a positive integer valued index. The basis 
functions chosen in EqU] are convenient when addressing 
biased conditions for the system. 

Note that left and right propagating electrons with the 
same s-index are not degenerate, but time reversal sym¬ 
metry is satisfied i.e. = E~^^ (simultaneous change 
of A and s). This symmetry reflects the fact that the SO 
interaction is not symmetric under space inversion but 
preserves time reversal symmetry (simultaneous change 
of A and s) so we retain Kramers degeneracy. The chiral¬ 
ity label C also reflects inversion asymmetry as changes 
in the chiral sign, at fixed A and s, change the energy. 


When the spin-orbit interaction is absent (a = 0), space 
and spin inversion symmetries are recovered (four fold 
degeneracy) once one combines fi = n/2N and h -|- 1 
labelled eigenvalues. 

Figure [T] shows a sequence of levels starting from the 
ground state (Kramers doublets) at two successive values 
of n along with their degenerate, at a = 0 partners with 
index -|-1, as a function of the SO strength. The ordering 
of the levels are indicated according to the spin orienta¬ 
tion and sense of the current of the chirality () = -1-1. 

To obtain a physical intuition on the nature of the 
wave functions (l^. , we derive explicitly some of the 

coefficients in the ansatz put forward. We explicitly do 
the A = -b, C = + case. Using equations dH) and (|4|) we 
find from the secular equation 


B 


+ _ 
-l-,s 


OJQ 

WSO 




( 6 ) 


where cos 6 = 1/yjl -b (wso/^o)^- In order to conform 
to a normalised spinor we choose Aljl ^ = cos(0/2), and 
thus ^ = sin(0/2), so we have 


tan 6 = 


i^SO 

Wo 


(7) 


This angle results from the existence of a SO “magnetic 
field” [illiH Hso = —a{k X z){2c/e) = {2ca/ahe)X\L\p, 
where \L\ is the angular momentum of the electron on 
the helix, and A gives the direction of the effective field 
depending on the rotation sense of the electron. 

The choice for the second eigenfunction is _ = 
— sin 0/2, leading to the two eigenspinors 


^+>+ _ 
^n,+ ^ 


/cos |e-*^/ 2 \ 
y sin|e®‘'’/^, y ’ 

/- sin|e-*^/2\ 

cos|e*‘^/2 J 


( 8 ) 


The corresponding eigenfunctions for A = — are 


_ -ihip 1 

/cos I 


n,+ — 6 1 

/ sin §e *‘^/2 J ’ 


^ g-inc, 1 

/- sin |e-*^/ 2 \ 

/ cos (e^v /2 J ■ 

(9) 


The angle 9 here tells about the inclination of the spinor 
with respect to the vertical z axisfl^ (see FiglT]) and 
n = n/2N with n an integer. When the SO coupling is 
very strong, 9 — > 7r/2, and the spin is in the plane, while 
for very small SO coupling, the spin is aligned with the 
z axis. 


III. SPIN FILTERING AND SUPPRESSED 
BACKSCATTERING 

The rotation and spin eigenvalue corresponding to the 
wave functions is depicted in Figl2]for both chiral labels. 
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FIG. 2. Electron propagation direction coupled to spin ori¬ 
entation and the corresponding energies for the two possi¬ 
ble chiralities of the helix. Each energy is doubly degenerate 
(Kramers doublet). If the electron propagates on a (" = -fl 
helix in the positive direction (as indicated), the lower en¬ 
ergy state corresponds to the spin up state, while the spin 
down state has to pay the energy gap A. 


selection of a propagation direction by applying a bias 
(as in / — y experiments) or by local population im¬ 
balance, effectively breaks time reversal symmetry by se¬ 
lecting one of the time reversed partners and generating a 
net spin collected on the other end of the chiral molecule. 
Note also that scattering in the helix is suppressed by this 
gap if the spin is not concurrently flipped by the scatter¬ 
ing event |22|. This is exactly the same mechanism for 
protected transport in chiral edge state in graohenef^. 

To further connect with the experimental results, we 
note that if we change the bias direction for the same 
chirality label, the preferred spin is the opposite spin di¬ 
rection by choosing the other partner of the Kramers 
doublet of lowest energy. This change of selected spin is 
proven in the experiments by the fact that spin injected 
back into the Ni magnet has to go where the DOS is 
higher, which is the opposite spin state [l^. 

The symmetry of the I — V curve in experiments is 
directly related to the fact that the same degenerate en¬ 
ergy corresponds to the forward and backward bias, the 
change only being which of the two partners in the dou¬ 
blet is selected. The barrier for electrons to be injected 
into the molecule depend on the workk function of the 
metal, and this is also very similar between Ni and Au 
(~ 5 eV), so this is also a source of symmetry. 


We have also evaluated the energies corresponding to the 
split doublets separated by the gap A. For fixed chirality, 
say C = +1, there is a doubly degenerate (changing s and 
A simultaneously) low energy configuration and a high 
energy configuration for each fi. 


E 


> 


fiujo 

2 


E 


< 


2 



( 10 ) 


As shown in FiglJl electrons propagating in the positive 2 
direction will have a lower energy if their spin is s = -1-1. 
The same energy corresponds to their Kramers partner 
which propagates in the opposite direction. In order to 
propagate a s = — 1 state in the positive direction, we 
need to pay an energy price A 


A = 


1 


n + - I 

2 ) 2ftwo’ 


( 11 ) 


for tuso/wo < 1 and Eqo = ^SO- On the other hand, 
if the energy Hloq < fiwso due to poor mobility of the 
electron on the helix (small ujq due to large effective mass) 
then the energy gap is directly related to the SO energy 

A = + 2^ (12) 


Thus, establishing a sense of electron propagation or bias 
will select a preferred spin direction (lowest energy). The 


IV. CHARGE AND SPIN CURRENTS 


Another way to pose spin selectivity is by evaluating 
the spin currents through the helix. In a coherent regime 
one can assess spin transport in the helical molecule by 
computing the expectation value of the velocity operator 
Jcharge = 'f'leu 4', where e is the electron charge and v 
is the velocity operator. In the presence of the spin-orbit 
interaction, the velocity operator is not simply p 1/m (a 
diagonal matrix) as there arises an additional anomalous 
velocity term. We start from the quantum mechanical 
definition of the velocity v = ^[H,r|. The azimuthal 
velocity component (/, is then 


—ihd^l2x2 — maacTp 


(13) 


The Hamiltonian takes a simple form when expressed in 
terms of the velocity: H = im*v^. This result is a man¬ 
ifestation of the possibility to write the SOC as a gauge 
fieldfij and thus a gauge invariant velocity as defined in 
EqUSl Note that all the manifestations of spin filtering 
have been observed in biased setups. In the free electron 
case, photoelectrons are emitted from the metallic sur¬ 
face in contact with one end of the molecules, and travel 
in a preferred direction. Also, in the localized regime, 
molecules are biased in a particular direction selecting 
the preferred sense of transport of the electrons. This 
bias, as argued before, breaks time reversal symmetry 
and chooses between the states depicted in FiglU We 
compute the charge currents for the biased system in the 
positive z direction 
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j+’+ 

'^n+l,+ 


hb 



m*y/a^ + 6 ^ 


(14) 


versus the opposite spin state in the same propagation 
direction 


hb 


A+ 

) 

2 

m*y/a? -I- 62 


(15) 


Thus, different spin states propagate at different ve¬ 
locities generating a net spin filtering effect. The net 
spin current vanishes for either SO zero or zero chirality 
(& = 0 ). ^ 

The difference is a net spin up current in the positive 
z direction 


J = J, 


+,-i- 


T+,+ 


_ p 

n+l,+ 


of magnitude from planar graphene to carbon nanotubes, 
just by bending the graphene sheets. In graphene, the in¬ 
trinsic SOC is second order in the atomic coupling {fieV 
second neighbour coupling) while it is linear in the atomic 
coupling for single walled nanotube (meV from direct 
Px{i)—Pz{j) coupling [3[ 2 ^ where i, j are nearest neigh¬ 
bours). The latter situation is also present in DNA if one 
regards the available orbitals to be on it bases, which 
are endowed with available Pz orbitals at such angle as 
to generate both V),p^ and Vppa couplings and linear in 
atomic SOC. The derivation is beyond the scope of the 
present paper, but the results show nontrivial features 
related to structural effects. 

We computed the effective coupling between DNA 
inter-base nearest neighbour Pz{i) orbitals, assuming one 
orbital per base, and a single helix. The same features 
should apply to the double helix. We found the following 
form for the SO coupling 


Vso oc 


2A(fyp - V^p) 

- e?p)^ ’ 


(18) 


hb 



(16) 


The longitudinal spin transport current is well defined 
and can be calculated as 




n+1,-1- 






2m*\/ 


62 


(cosd — 1) . 


(17) 


The longitudinal spin current depends on the spin-orbit 
coupling through 9. The spin filtered in the 2 ; direction 
disappears when the SO coupling is zero (0 = 0). Again, 
the pitch dictates the strength of the vertical spin cur¬ 
rent and both the pitch and the SOC must be present. 
Note that the factor (1 — cos0) = is the non- 

adiabatic Aharonov-Anandan phase found by Frustaglia 
and Richter in a detailed analyses of conductance through 
SO coupled rings [l^. This phase offers a new insight into 
spin filtering of chiral molecules since in the strong SO 
limit it is related to the Berry phase of the spin and for 
the general case it is controlled by non adiabatic spin 
precession. 


V. STRENGTH OF THE SOC AND ROOM 
TEMPERATURE SPIN SELECTIVITY 

Recent experimental findings demonstrating spin selec¬ 
tivity have shown that the gap A in EqlT^l is much larger 
that expected from even the purely atomic contribution. 
We have attempted to explain this anomalous size of the 
SO gap through a tight-binding approach[T^ [l^, It 
is well known that the SOC can be enhanced three orders 


The parameter A is the bare atomic SOC, Vpp and Vpp 
are the Slater-Koster p orbital couplings perpendicular 
and parallel (respectively) to the line joining the carbon 
atoms of the model, ejp — i® the difference in en¬ 
ergy between the bonded p^ y in the plane of the base 
and the Pz orbital perpendicular to the plane. Note that 
the energy contrast in the denominator may be small, 
potentially increasing the SOC by a large amount. This 
is especially true when interbase coupling is weak. An¬ 
other way to express the denominator using lowest order 
perturbation theory is by renormalizing the in plane cou¬ 
pling by the s -p overlap; - ejp - (KpV(e 2 p - £ 2 *), 
where e^i are the bare atomic values and V^p is an in¬ 
terbase coupling that needs to be computed from a first 
principles study. 

The spin selectivity demonstrated by experiments is 
a room temperature phenomenon. Thus any coherent 
mechanism must be limited by decoherence lengths [T^. 
M- The gap for degradation of spin selectivity found 
in experiments is very high, 0.5 eV (refl3), compared 
to thermal effects at room temperature 25 meV. This 
gap prevents elastic backscattering and exponentially re¬ 
duces inelastic scattering with the same spin. A spin cou¬ 
pled scattering mechanism, nonetheless, could degrade 
the spin selectivity rapidly because of the existence of 
the Kramers doublets as a backscattering channel (see 
Figlll). 

If one proposes a decoherence length operating at 
room temperature, one can also suggest a mechanism 
which preserves spin selectivity analogous to that pro¬ 
posed in ref. 0 and 0, where relaxation of phase co¬ 
herence occurs within a few nanometers while spin se¬ 
lectivity is preserved. This mechanism entails an exp o- 
nential decay of the transmission with a decay rate[^ of 
P = -(l/dcoh)ln[teff/(Aa-(A_F-eU/2))] where Aa is the 
energy of the spin preffered channel, teff is the hopping in¬ 
tegral between sites separated by a coherence length dohi 




















6 


Ep is the Fermi energy and V the applied bias. Then in¬ 
creasing the number of dcoh by having longer molecules 
will make the input current decrease exponentially, so 
and increased bias is needed to achieve the same out¬ 
put current. This mechanism might be the source of the 
increased barrier for spin selectivity as the molecule is 
elongated. 

Again, there is no mechanism for the coupling to 
backscattering channels with the same spin (See Figl2|) 
so that there will be degradation of spin current but not 
of spin polarization. 

VI. SUMMARY AND CONCLUSIONS 

We have examined a model for chiral spin selectivity 
on a spin orbit active helix with N turns. The origin 
of the SO electric field comes from the atomic cores of 
the carbon atoms which provide, through the Pz orbitals, 
a narrow band for transport. The resulting channels of 
the model helix are Kramers doublets involving opposite 
propagating and opposite spin projections at the same 
energy. An applied bias or otherwise preferred transport 
direction will then select a spin and effectively break time 
reversal symmetry (choosing one in a doublet pair), and 
transferring a particular spin. The resulting spin cur¬ 
rent is the spin selectivity mechanism. The spectrum of 
the model also insures suppressed backscattering by an 
energy gap controlled by the SO energy and the inter- 
base coupling of the Pz orbitals, as suggested by a tight 
binding calculation. A form for the coupling is proposed, 
that could enhance the SOC in lowest order perturbation 
theory. 

There are a few predictions of the theory that could be 
tested experimentally: 1) The spin selected, after chang¬ 
ing the chirality of the transport structure, seems to be 
the same although the state is different. 2) There should 
be changes in the strength of the SOC (as measured 
through the experimentally measured gap), by DNA de¬ 
formations that alter the inter-base hopping integrals. 
These could be induced by changing torsion (tugging on 
the molecule). The SOC, according to the tight-binding 


model, should increase as long as the molecule is not dis¬ 
rupted. 3) Current decay at fixed bias when changing 
the molecule length. 4) Depending on the strength of 
the SOC and the length of the molecule, the number of 
same-spin channels can increase making the selectivity 
more robust. 5) The suppressed backscattering mecha¬ 
nism could be tested by changing DNA bases introducing 
random defects that do no couple to spin. The spin se¬ 
lectivity should be robust against such inhomogeneities. 

The model here can be made more quantitative in 
many directions. The derivation of the Hamiltonian di¬ 
rectly from the tight-binding descriptions is desirable so 
that one can properly account for the geometry of the or¬ 
bitals participating in transport. We have solved a vari¬ 
ation of the problem more akin to the geometry of the Pz 
orbitals in carbon nanotubes, where they rotate on the 
outside of a cylinder. Although the resulting Hamilto¬ 
nian is different in detail the same physics described here 
follows. This is expected from the symmetry arguments 
that result in Kramers doublets separated by a gap. Con¬ 
templating the double helix structure of DNA, could also 
bring about new interference effects that might enhance 
or reduce spin selectivity [l^ and bring the model quan¬ 
titatively closer to the experiments. 

Finally, incorporating the metallic contacts to the chi¬ 
ral molecule through the thiol groups, should introduce 
level broadening to the spectrum of the molecule and de¬ 
termine the escape rate of electrons and possible charging 
phenomena. 
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